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Abstract — We present a general dyadic vector circuit formalism, applicable for 
uniaxial magneto-dielectric slabs, with strong spatial dispersion explicitly taken into 
account. This formalism extends the vector circuit theory, previously introduced 
only for isotropic and chiral slabs. Here we assume that the problem geometry 
imposes strong spatial dispersion only in the plane, parallel to the slab interfaces. 
The difference arising from taking into account spatial dispersion along the 
normal to the interface is briefly discussed. We derive general dyadic impedance 
and admittance matrices, and calculate corresponding transmission and reflection 
coefficients for arbitrary plane wave incidence. As a practical example, we consider a 
metamaterial slab built of conducting wires and split-ring resonators, and show that 
neglecting spatial dispersion and uniaxial nature in this structure leads to dramatic 
errors in calculation of transmission characteristics. 
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1. INTRODUCTION 

It is well known that many problems dealing with reflections from multilayered 
media can be solved using the transmission line analogy when the eigen-polarizations 
are studied separately (see e.g. p^). In this case, the amplitudes of tangential 
electric and magnetic fields are treated as equivalent scalar voltages and currents 
in the equivalent transmission line section. In order to account for an arbitrary 
polarization, Lindell and Alanen introduced a vector transmission-line analogy |2j, 
where vector tangential electric and magnetic fields serve as equivalent voltage 
and current quantities. Later on, the vector transmission-line analogy was further 
extended for isotropic and chiral slabs into a vector circuit formalism, with the slabs 
represented as two-port circuits with equivalent impedances and admittances [31 13 • 
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This vector circuit theory has been successfully applied to study plane wave reflection 
from chiral slabs [Sj, and extended to uniaxial multilayer structures 

Recent emergence of metamaterials and the subsequent growth of research 
interest to their properties, revitalized the importance of analytical methods for 
studying artificial media, in particular, for proper calculation of the reflection 
and transmission properties. The most prominent examples of metamaterials [Zj 
involve split-ring resonators |Hj interlayered with a wire structure, being organized 
in an essentially anisotropic manner (e.g., like shown in Fig. ^,). Effective 
permittivity (permeability) tensors in such metamaterials correspond to those of 
uniaxial dielectric (magnetic) crystals, and the principal values of both tensors differ 
from unity along one direction only. Moreover, the presence of wire medium imposes 
significant spatial dispersion for all waves with an electric field component along the 
wires 

These peculiarities drive the properties of such metamaterials far apart from 
what can be described in terms of the isotropic vector circuit theory, and raise a 
clear demand for an appropriate generalization. The corresponding theory extension 
is the main objective of this paper. 



2. TRANSMISSION MATRIX 



Let us consider a spatially dispersive slab having thickness d and characterized by 
the following material parameter dyadics (Fig. ^>) 

e = e t (kt)7 t + e n nn, p = ^7 t + n n na, (1) 

where n is the unit normal vector for the slab and I. = I — nn is the transversal unit 
dyadic. We consider plane wave excitation and move from the physical space to the 
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Figure 1. a) A TM-polarized plane wave incident on a slab implemented using an 
array of wires and split-ring resonators, b) Macroscopic representation of a uniaxial 
magneto-dielectric slab. Subscripts + and — denote the fields at the left and right 
sides of the slab, respectively. 
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Fourier space by a transformation Vt — > — jkt. For plane waves, kt simply stands 
for the transversal propagation factor. Physically we could as well assume that the 
slab is illuminated by a source located electrically far away from the slab. Notation 
e(fct) stresses the spatially dispersive nature of the slab in the tangential plane, 
and indicates the dependence of the permittivity component from the tangential 
propagation factor. 

Starting from the Maxwell equations the following set of equations can be 
derived for the tangential field components 

d 1 

— n xE ( = -juntHt + k t x k t x Ht, (2) 

dz JUJ€ n 

d 1 

— n x H 4 = juje t (k t )E t - k t xk,x E f . (3) 

dz JUH n 

Next, we integrate eqs. (J5J) and Q over z from to d: 

n x E+ , — n x E* ~ 1 ~ 

±— = -jufitHt + k t xk(X H t , (4) 

d jioe n 

n x H> , — n x EL N ^ 1 ~ 

^— = jue t (k t )Et - k t x kt x E t . (5) 

d JLOfln 

Above E t+ ,H t+ refer to the fields at the left side of the slab, and E f _,Hj_ refer to 
the fields at the right side of the slab, Fig. ^i. The averaged fields in eqs. (@J) and 
© are defined as 

1 f d „ . - 1 



B t = -J E t dz, H t = -J U t dz. (6) 
After mathematical manipulation Q and (JSJ) transform into 

E t _ - E t+ = -ju/itdA • n x H 4 , (7) 

Ht_ - Ht + = jioe t (k t )dB • n X E t . (8) 
Above dyadics A and B are defined as 

= = = ktkt = io 2 ^n ~ kt kf kf n x k f n x kt 

* io 2 /j, t e n oo 2 {i t en k% k 2 

= _= _ k f k t _ uj 2 fi n e t {k t ) - fc| kt k t + n x k t n x k t ^ 

1 uj 2 fi n e t (k t ) to 2 fi n et(k t ) kj k\ 

The general solution for the transverse electric field inside the slab reads (now the 
interest lies only on wave propagation in z-direction) 

Et(z) = A ■ e~^ z + B ■ e^ z , (11) 

where A and B are constant vectors and the z-component of propagation factor is 
different for TM and TE polarizations: 



kt kt n x kt n x k 



(3 = (3(k t ) = (iju—r-T + /?te 72" 

fc t fe t 



et(kt) , 2 2 ktkt / // t f o ,, x 2 . n x kt n x kt 



(^*e„ _ fc 2 } ^t + i^^ et{h) _ k 2 t) t- ^ (12) 
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Constant vectors A and B are determined from the boundary conditions 

E(0) = E t _, E(d) = E t+ , (13) 
and the following expressions can be derived 

A = {e~^ d - e^ d Y l ■ (E t + - E t _ • e^ d ), (14) 

B = -(e-^ d - J?*)- 1 ■ (E t+ - Et_ • e~^ d ). (15) 

After integrating 1)11(1 over z from to d we get for the averaged electric field (note 
that all the dyadics are commutative) 

E t = —(E t+ + E t _) • /I" 1 ■ (e~^ d - e^ d y l ■ (2% - e~^ d - e^ d ) 
3d 

= _ L tail (¥) 
= (E t+ + E t _) • /3 y^L. (16) 



Similarly for the magnetic field 

tan 



H t = (H t+ +H t _)-/3 (17) 

a 

Inserting ()16|) and ()17() into (J7J) and @ leads to the following result 

E t _ - E t+ = -jujfitdl ■ n x (H t+ + H t ) • T 1 ■ 7, (18) 

Ht_ - H t+ = juje t (k t )dTi • n x (E t+ + E t _) ■ /T* ■ 7- (19) 

Above we have denoted / = tan((3d/2)/d. Fields at the upper side of the slab can be 
expressed with the help of the fields at the lower side of the slab and the following 
result is obtained after mathematical manipulation: 

E t+ = cos 0d) ■ E t _ + jufjLt sin (fid) ■ /T • I • (n x H t „ ), (20) 

n x H i+ = cos (j$d) • (n x H t _ ) + j— sin 0d) ■ (jf 1 • IV 1 • E t _ . (21) 

tout 

Writing (|2U|) and 1)21(1 into a matrix form we identify the dyadic transmission matrix 
for the slab, Fig. 



E t+ \ = ( a u «uj V f E *- 
n x H t+ J V a 2i a 22 / V n x H * 



(22) 



where the transmission components read 



kf kf n ,,nxkf nxkj 



a n = a 22 = cos(P TM d)—^Y- + cos(/3 TE d) ^ , (23) 

t 



V t - ' k 



= .uje t (k t ) . „k(k( ./3 TE . n x k 4 n x k t 

a 2i =J— S sm{P TM d)— — + j sm(/3 TE d) , (24) 

Ptm w//f K t 

wet fct fc£ Pte kf 
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We immediately notice that if the slab is local and isotropic, coefficients (|23|) - (|25|) 
reduce to those obtained earlier for an isotropic slab 

The exact boundary condition for a slab on a metal ground plane follows directly 
from (H3) with E 4 _ = 0: 

E t+ = I(t t )-nxH (+1 
where the impedance operator reads 



Z[kt) =3 7rT tan (PTMa)^2- + J -5— tan(/^ TE d) . 



(26) 
(27) 



3. IMPEDANCE AND ADMITTANCE MATRICES 

From (|22|) it is straightforward to derive the impedance and admittance matrices for 
the slab 

' B t+\ = ( Z_n 1u\ ( n x H t+ 



n x H t+ 
n x H 4 _ 



y y 

£.11 _j_12 
^21 ^22 



Et+ 
E t _ 



(28) 
(29) 



where the dyadic impedances and admittances depend on the transmission 
components in the following way 



Z n — a u ■ a 21 , Z 12 — - 
Z 21 = a 21 , Z 22 



= — a 21 ■ a 22 , 



y il — a 22 ■ a l2 , Y 12 — a 21 a 22 ■ a l2 ■ 



Y 21 — a 12 , Y 



22 



-a 12 -a u . 



(30) 
(31) 

(32) 
(33) 
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Figure 2. Different representations for the slab: a) A two-port transmission line 
section, b) Vector T-circuit. c) Vector ri-circuit. The orientation of the unit vector 
n is in all the cases the same as depicted in a). 
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The corresponding T and l~l-circuit representations are presented in Fig. and 
Fig. [2fc, respectively. 

4. REFLECTION AND TRANSMISSION DYADICS 

Introducing dyadic reflection and transmission coefficients R and T, eq. (|22[) can be 
written as two equations in the following form: 

Ch + 1) • Ej nc = 5 U • f ■ Ef c + i 12 • 1 o 1 • f • E» c , (34) 

f „ 1 • (f t - 1) • e» ic = * 21 . f ■ Er + s 22 • i o 1 • ? . e» ic , (35) 

where E mc denotes the incoming electric field and Z is the free space impedance 
dyadic (seen by the tangential fields) 

= k«kt 10 nxk |n xkj 

Z ° = %cos0-^- + — . (36) 

From eqs. (|34)1 and (|35|) we can readily solve the transmission and reflection dyadics: 



T = 2(a n + a 22 + a 12 - Z + Z -a 21 ) 1 • I t , (37) 
i? = (5 U +ai 2 • Z" 1 ) -T-7 t . (38) 
5. PRACTICALLY REALIZABLE SLABS 

A typical example which fits into the general model presented above, is a slab of 
metamaterial, implemented as an array of conducting wires and split rings resonators 
(WM-SRR structure), Fig^i. The wires are assumed to be infinitely long in the 
y-direction. Moreover, the number of wires in the x-direction, and the number of 
split-ring resonators both in the x- and y-directions is assumed to be infinite. 
For the case shown in Fig. the non-local permittivity dyadic reads [2] 

e = e x u x u x + e y u y Uy + e z u z u z , (39) 

e x = e z = e h , e y = e h (l - J , (40) 

where is the permittivity of the host matrix, k p is the plasma wave number, k 
is the wave number of the host medium, and k y is the y-component of the wave 
vector inside the lattice. Generalization of the dyadic for a case when the wires are 
perodically arranged along two directions (double-wire medium with non-connected 
wires ^0]) is straightforward. 

A commonly accepted permeability model as an effective medium description of 
dense (in terms of the wavelength) arrays of split-ring resonators (SRRs) and other 
similar structures reads (see e.g. [TT1 IT21 IT3*] ) 

]I = fJ, X U X U X + fJ,yUyUy + fJ, Z U Z U Z , (41) 

1 + ~~2 2, t ) ' A*w = = Mh, (42) 
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where fi^ is the permeability of the host medium, A is the amplitude factor 
(0 < A < 1), ujq is the undamped angular frequency of the zeroth pole pair (the 
resonant frequency of the array), and V is the loss factor. These parameters can be 
theoretically estimated for any particular case When the SRRs are positioned 
to the locations where the quasi-static magnetic field produced by the wires is 
zero (to the symmetry planes), there is no near field coupling between the two 
fractions |14j and the whole metamaterial can be characterized by the permittivity 
and permeability in the form above. 

The derivation presented above holds for the WM-SRR structure only when 
the wires are parallel to the slab interfaces and the tangential propagation factor is 
restricted to the interval < kt < k. In this case a plane wave incident on the slab 
excites only TM and TE wave. If the slab is excited by a source, located close to 
the surface, or the source is inside the slab, an additional TEM wave will be excited 
by the source, [§]■ The same happens for an incident plane wave, if kt = k. Note 
that, when the wires are perpendicular to the interface different approach is needed. 
Recently, authors of considered a slab with wires perpendicular to the interfaces 
and presented another method to calculate the transmission coefficient for such a 
slab. 

6. SPECIFIC EXAMPLES 

Here we present some calculated results for plane wave transmission through the 
WM-SRR slab in Fig. ^p. We compare the transmission coefficient calculated when 
the slab is assumed to be uniaxial and non-local with the transmission coefficient 
obtained when: (i) the slab is assumed to be uniaxial but local, or (ii) the slab 
is assumed to be isotropic and local. Local model for the permittivity means that 
k y = in (|40j). Next, we study the angular dependence of the transmission coefficient 
at certain frequencies when the effective refractive index of the slab is close to zero. 
For simplicity, only an example of TM polarization is considered. 

For this analysis we assume the following parameters: slab thickness d = 150 
mm, eh = eo> an d k p = 104.7 m _1 (the corresponding plasma frequency is / p = 5 
GHz); = /^OjA = 0.4, u^o = 2-7T • 2.5 GHz, T = ujq/50. The permittivity and 
permeability dyadics are calculated in accordance with equations ()39j) - ()4U|) and (|41[) - 
(|4"2|) . respectively. 

Fig. |31 compares the exact transmission coefficient to the transmission coefficient 
calculated when the slab is assumed to be local [case (i)]. Only for the normal 
incidence (not plotted) the results are the same. At small incidence angles (e.g. 9 = 
7r/6) there is a transmission maximum around 3 GHz. This maximum corresponds 
to the frequency range where both Re{e} and He{fi} are negative and relatively 
close to unity in magnitude. In this situation a planar slab can bring a point source 
to a focus without spherical aberration ^H]- In a certain frequency above 3 GHz 
Re{/j,} becomes positive while Re{e} remains negative, leading to a stop-band. At 
frequency f = f p permittivity becomes positive and waves can propagate through 
the slab. Note that when the permittivity is assumed to be non-local the position of 
the pass-band edge is predicted at remarkably higher frequencies compared to the 
local model. 

For larger incidence angles (e.g. 9 = tt/3) the maxima around 3 GHz disappear 
and the transmission coefficient obeys an increasing behavior starting at /o (the 
resonant frequency for fj, t ). Indeed, eq. (|T2*)) for (3 JU shows that the term Rej/c^^f — 
sin 2 #)} is negative at frequencies / > /o- Thus, a pass-band will appear in the 
frequency interval fo<f<f P because in this interval also Re{e} is negative. 

Fig. ^compares the exact transmission coefficient to the transmission coefficient 
calculated when the slab is assumed to be local and isotropic [case (ii)]. Clearly the 
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F [GHz] 

Figure 3. Transmission coefficient (absolute value) as a function of frequency at 
different incidence angles. Exact transmission coefficient (thick lines) is compared 
to the transmission coefficient calculated using local permittivity model (thin lines). 
Solid lines: 6 = ir/6, dashed lines: 9 = ir/3. 

assumption that the slab is isotropic and local leads to severe errors except for the 
normal incidence. 

For an additional illustration we study the angular dependence of the 
transmission coefficient at certain frequencies when the effective refractive index 
n of the slab is close to zero while Re{/^} > A Re{fit} ^ A*o> an d Re{e t } > A 
Reje^} ^ eo- Comparison is made between the results given by the non-local and 
local models (in both models the slab is assumed to be uniaxial). An interest to this 
problem arises from recent suggestions for microwave applications benefiting from 
slabs having low value for both Re{/i} and Re{e} [T7] . 




F [GHz] 

Figure 4. Transmission coefficient (absolute value) as a function of frequency at 
different incidence angles. Exact transmission coefficient (thick lines) is compared to 
the transmission coefficient calculated using local and isotropic permittivity model 
(thin lines). Solid lines: 6 = tt/6, dashed lines: 8 = tt/3. 
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Figure 5. Transmission coefficient (absolute value) as a function of the incidence 
angle at certain frequencies. Exact transmission coefficient (thick lines) is compared 
to the transmission coefficient calculated using local permittivity model (thin lines). 
Solid lines: / = 1.05 x f p , dashed lines / = 1.20 x / p . 

In practise the condition n « 0, Re{/i(} > A Reje^ > occurs at frequencies 
slightly above / p . We have to bear in mind, however, the physical limitation of the 
permeability model (|42|): the model is valid at low frequencies and at frequencies 
relatively close to the magnetic resonance. Here we assume that in the vicinity of 
fp the permeability model is still valid. The transmission coefficients are depicted 
in Fig. H3 It shows that transmission maxima, covering a wide range of angles from 
50 to 80 degrees appear when the non-local model is used. These maxima are not 
predicted when the local permittivity model is used. This behavior can be explained 
by considering the z-component of the propagation factor, which in this case can be 
written in the following form: 

In order for a wave to propagate through the slab, both terms inside the parentheses 
[the right side of eq. ([43))] must be simultaneously positive or negative. For small 
incidence angles both are positive and for large angles both are negative. In a certain 
range of angles, however, these terms are of opposite sign, leading to a stop-band. 

Effectively condition n ~ can also be achieved in photonic crystals in the 
vicinity of the stop-band edge, e.g. ^3 EE HH- This feature is reported to be 
important for practical applications (e.g. |U). Accordingly, we can apply the 
developed method to study the transmission characteristics of the slab only in the 
presence of wires. Fig. 03 shows the calculated results (note the range of incidence 
angles). The transmission maxima seen at certain angles correspond to thickness 
resonances of the slab. The results indicate that the slab can be utilized as an 
effective angular filter at microwave frequencies. Note that when the permittivity 
is assumed to be local, some maxima are also seen in the transmission coefficient, 
however, the location of these maxima is incorrectly predicted. 
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Figure 6. Transmission coefficient (absolute value) as a function of the incidence 
angle at certain frequencies when the slab is assumed to consist only of wires. Exact 
transmission coefficient (thick lines) is compared to the transmission coefficient 
calculated using local permittivity model (thin lines). Solid lines: / = 1.05 x / p , 
dashed lines / = 1.20 x f p . 

7. CONCLUSIONS 

In this paper we have formulated a vector circuit representation for spatially 
dispersive uniaxial magneto-dielectric slabs. A dyadic transmission matrix and the 
corresponding impedance and admittance matrices have been derived. The results 
take into account spatial dispersion along the planes parallel to the slab interfaces. 

The presented results allow the exact calculation of the transmission and 
reflection coefficient for a plane wave with arbitrary incidence angles. This model 
is applicable, for example, to a typical metamaterial implemented as a lattice of 
conducting wires and split-ring resonators. It has been shown that for accurate 
transmission analysis the uniaxial nature of such a slab, and the spatial dispersion 
in the wire media must be taken into account. The calculated results also indicate 
the feasibility of the slab to operate as an effective angular filter at microwave 
frequencies. 
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